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The dynamics of a two-state system whose energies undergo a real crossing at some instant of time 
is studied. At this instant, both the coupling and the detuning vanish simultaneously, which leads to 
an exact degeneracy of the eigenenergies of the system. It is found that the dynamics of the system 
is primarily determined by the manner in which the degeneracy occurs. This interesting behavior 
is reminiscent of a symmetry breaking process, since the totally symmetric situation occurring at 
the crossing is significantly altered by infinitesimal quantities, which remove the degeneracy, with 
very important dynamical implications from there on. A very simple analytical formula is derived, 
which is found to describe the population changes very accurately. 
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I. INTRODUCTION 

In the realm of adiabatic evolution [l|, the dynamics 
of a two-state system whose energies cross at some point 
plays an important role. The exactly solvable model 
of Landau-Zener-Stueckelberg-Majorana (LZSM) is 
the most famous example of such class of problems. Its 
widespread applications range from quantum optics and 
atomic physics to chemistry and biophysics. In fact, this 
model is of interest in several physical systems [^, such 
as Josephson junctions 0, cold atoms in optical lat¬ 
tices @ and spinorial Bose-Einstein condensates M- In 
order to better consider realistic situations, the implica¬ 
tions of having a finite duration or a nonlinear 

crossing [I^ in the experiments have also been studied. 
Generalizations of the original model to the cases wherein 
more than two states are involved have been considered 
[30, as well as level crossings occurring in systems 
described by nonlinear equations El- The relevance of 
the LZSM model in statistical physics, especially in con¬ 
nection with phase transitions in the Ising model, has 
been discussed (30- Moreover, recently, a geometric 
Landau-Zener-Stiickelberg interferometer has been real¬ 
ized in single trapped ions by Zhang et al [3| on the 
basis of a proposal by Gasparinetti et al [3 | , and Landau- 
Zener transitions in frozen pairs of Rydberg atoms have 
been observed by Saquet et al 

In the original LZSM model, the bare energy levels of a 
two-state system cross at some instant of time, but, due 
to the presence of a (constant) interaction, the degener¬ 
acy is removed in the adiabatic eigenenergies. Then the 
adiabatic following of an instantaneous eigenstate of the 
Hamiltonian leads to a perfect transition from one bare 
state to the other. Beyond the perfect adiabaticity, the 
correction to the dynamics is given by the formula found 
by Landau, Zener, Stiickelberg and Majorana which ex¬ 
presses the transition probability between the eigenstates 
(non-adiabatic transitions) of the Hamiltonian , and 
which is equal to the no-transition probability between 
the bare states. 

In 1990 Fishman et al have introduced a variant of 


the original LZSM model where neither the bare ener¬ 
gies nor the dressed ones really cross 0], showing that, 
even in the presence of the same dressed energies, the 
circumstance that the bare energies do not cross as in 
the standard model determines significatively different 
results. One of the most interesting experimental re¬ 
alizations of this modified LZSM model has been pro¬ 
vided by Bouwmeester et al , who have observed such 
modified Landau-Zener dynamics (they call the relevant 
mathematical description the ‘hidden-crossing model’) in 
classical optical systems. By the way, it is worth men¬ 
tioning that with their experiment Bouwmeester et al 
have demonstrated that the validity of the LZSM model 
(and related ones) is not limited to the realm of quantum 
mechanics. 

The variant of the LZSM model complementary to the 
hidden-crossing model corresponds to a physical situa¬ 
tion where both bare and dressed energies cross at the 
same time, which has never been studied. It is then of 
interest to analyze the adiabatic evolution of a two-state 
system in the special case where a real degeneracy occurs 
at the crossing of the bare energies. In other words, we 
want to study the case where at the same time when the 
bare energies cross, the interaction, assumed non van¬ 
ishing in the rest of the process, vanishes too. At this 
instant the adiabatic (i.e., dressed) energies, just like the 
diabatic (i.e., bare) energies, become degenerate. In this 
case the LZSM formula is not applicable, and in order 
to estimate the transition probability we have to scruti¬ 
nize the intricate behavior of the eigenenergies and the 
nonadiabatic coupling in the vicinity of this real crossing. 

The paper is organized as follows. In the next sec¬ 
tion, we introduce the model Hamiltonian we are going 
to analyze. In section Hill we discuss the condition for the 
validity of the adiabatic approximation, firstly in general 
and secondly focusing on the instant of time when the 
level crossing occurs. In section HVl we report on an anal¬ 
ysis that allows to predict the efficiency of population 
transfer on the basis of the peculiarities of the Hamilto¬ 
nian. Finally, in section |V] we summarize the results and 
provide some conclusive remarks. 
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II. THE MODEL 


Consider a two-state system whose dynamics is gov¬ 
erned by a Hamiltonian having the following form {h = 
1): 


„ A n 

H = — a, H-(T 

2 2 



-A n\ 
n A )- 


( 1 ) 


This Hamiltonian has eigenvalues 


E±=u]±= ± —Vl + 0? , 

(2) 

and eigenstates 

\il)+) = cos6l 1) -1- sinS f) , 

(3a) 

^_) = sin0 ),) — cosO t) , 

(3b) 

where 

tan 9 = a+ \/l + a'^, 

(4a) 

A 

“ H ■ 

(4b) 


In the LZSM model, the interaction is constant {Ct = 0) 
and the bare energies change linearly (A = nt). If the 
change is very slow (k small enough) then the evolution is 
adiabatic. Moreover, when t varies from —oo to -l-oo, the 
dressed states change from one bare state to the other; 
for example, |r/>+(—oo)) = H) becomes |'0+(oo)) = |t)). 
Therefore, if the system is prepared in one of the bare 
states, and the system follows the corresponding dressed 
state, then in the end of the process each bare state is 
perfectly mapped to the other bare state: population 
inversion has occurred. Perfect adiabatic evolution does 
not exist (except for constant Hamiltonian or for A = 0), 
and the LZSM formula gives the deviation from perfect 
population inversion. 

In the following, we will instead study the adiabatic 
evolution when H and A are such that A —?► 0 and —>■ 0, 
when t —>• 0, so that a real level crossing occurs at the 
time instant t = 0. Of course, there are several ways the 
two quantities can approach zero, and depending on this, 
the evolution of the system can be very different. 

In order to maintain the possibility of having a com¬ 
plete population transfer, we will always assume that a 
varies from —oo (i.e., a huge negative value) to -l-oo (a 
huge positive value) or vice versa. Indeed, this will make 
each of the dressed states change from being exactly (i.e., 
essentially) equal to one bare state to being exactly (i.e., 
essentially) equal to the other bare state along the exper¬ 
iment. 


III. ADIABATIC EVOLUTION 

After introducing the unitary operator describing the 
passage to the basis of the instantaneous eigenstates of 


the Hamiltonian, 

Ta = |^/’+(-oo)) (V’-h(t)l + IV’-(-oo)) (^-(i)l 

= ||)(^+(t)| + |t)(^-(t)l , (5) 

we easily find that in this basis the evolution is generated 
by the following operator: 

- iTATi+ taoti = ). (6) 

It is worth noting that for 0 = 0 one would seemingly 
have perfect adiabatic evolution (in fact, each of the 
two instantaneous eigenstates would just acquire a phase, 
without undergoing any transitions). But, if 0 = 0 then 
there is no change of the eigenstates of the Hamiltonian, 
and then there is no real adiabatic evolution. In fact, 
the adiabatic following of the eigenstates of the Hamil¬ 
tonian governing the system is intrinsically the result of 
some approximations and obtaining an exact adiabatic 
following of the eigenstates of the system Hamiltonian is 
a chimera. There is, seemingly, a violation to such an as¬ 
sertion given by the shortcuts to adiabaticity [25I - [T7I| . In 
fact, such techniques allow to realize the exact adiabatic 
following of the eigenstates of a given Hamiltonian. But 
the crucial point of such shortcuts is that the exact fol¬ 
lowing is obtained through an additional interaction, so 
that, in the end, the adiabatic evolution does not coincide 
with the adiabatic following of the eigenstates of the sys¬ 
tem Hamiltonian. In other words, in order to realize the 
adiabatic following of the eigenstates of a given Hamilto¬ 
nian H{t) it is necessary to let the system evolve under 
the action of H{t) + Hsc{t), where Hsc{t) is a suitable 
additional term that allows to realize the shortcut and 
which does not commute with H{t). In this way, there is 
no following of the eigenstates of the total Hamiltonian 
of the system H{t) + Hsc{t). 


A. General condition for adiabaticity 


In order to have an adiabatic evolution the terms re¬ 
sponsible for transitions between the eigenstates of the 
Hamiltonian must be negligible, which in our case means: 
r] = |0|/|w+ — w_| <C 1. Now, since 


Vl + 


I -|- (a -I- vT+T?)^ vT+T? ’ 

one finds: 

V = 

a + \/l + 


(7) 


1 -I- (a -I- Vl -I- a^Y 


H(I-ha2) 


, ( 8 ) 


which can also be cast in the following form (consider 
that a -I- vT+7? = l/(—a -I- vT+l?) and then 
2\J\ -I- o? = tan0 -I- 1/ tan0): 

T] = I sin(26l)|^ X I d/H| . 


(9) 
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Putting rj in this form, it is clear that a sufficient con¬ 
dition to guarantee an adiabatic evolution is | d/ff | 1. 
Nevertheless, in some cases, a finite or even a divergent 
value of I a/Vl\ can be compensated by the smallness of 
sin(2d). 

Let us come back to Eq. m and consider what the 
condition for adiabaticity becomes in the cases |q;| <C 1, 
a ^ 1 and |a| ^ 1. First of all, consider that both for 
|a| 1 and a ~ 1 one can say that rj ~ |d/n| = |(f7A — 

nA)/n^|, while for |a| ^ 1 one has ij ~ |d/(a^fl)| = 
|(f2A — flA)/A^|. Therefore, the conditions to have an 
adiabatic evolution become: 

{ |(f7A-ffA)/ff3|«i, (10a) 

|a|»l ^ |(f7A-f2A)/A3| < 1. (10b) 

Note that when |a| 1, for a > 0 one has 9 —>■ 7r/2, 

while for a < 0 it turns out 0 —>■ 0, and then sin(20) —>• 0 
in both cases. This is just the case mentioned above 
where smallness of sin(2d) can compensate higher values 
of |d/n|, and this is the reason why the condition for 
having an adiabatic evolution assumes a different form. 

B. Adiabaticity at the level crossing 

It is of interest to give some general indication of the 
validity of the adiabatic approximation close to the level 
crossing. In order to develop such an analysis, let us 
consider the behavior of A(t) and 11 (t) close to the level 
crossing occurring at t = 0, assuming that the leading 
terms in their Taylor expansions versus t at the crossing 
{t = 0) are 

A ~ , (11a) 

n~Bt\ (11b) 

with a,b > 0. Then, in the case a <C 1 (a > &), one has: 


IV. NONADIABATIC TRANSITIONS 

In this section we give some general dynamical proper¬ 
ties that allow to forecast the efficiency of the population 
transfer on the basis of an analysis of A(t) and at 
the level crossing. 

A. Transitions at the level crossing 

Let us first of all assume that the evolution of the sys¬ 
tem is essentially adiabatic everywhere, with the only 
possible exception of the very neighborhood of t = 0. 
The total evolution operator (that is, the evolution op¬ 
erator in the Schrodinger picture) can be written as: 

Ui+oo,-oo) = t/(-hoo,0+)C/(0+,0-)17(0“,-oo) 

« t/A(+oo,0+)C/(0+,0-)17A(0-,-oo), 

(14) 

where we have introduced the adiabatic evolution oper¬ 
ator: 

UA{t2,ti) = ^ |?/'fc(i2)) ('!/'fe(<i)| , (15) 

k=± 

with 

= f (wk{s) - i (^i{jk{s)\'tpk{s)\'] ds . (16) 

Jti 

Now, since around t = 0 the Hamiltonian H does not 
exhibit singularities (far from this, it is infinitesimal), one 
has t/(0+, 0“) ~ I, and then: 

17(-l-oo, —oo) « t/A(+oo, 0''‘)17a(0“, —oo), (17) 

from which: 

17(-l-oo, —oo) = EE Tki |V’fc(+oo)) (18) 

A:—± l=zL 


V 


A(b — a) 




x\t\ 


i-2b-l 


and to guarantee ry <C 1 one needs a > 2b+1. 
In the case a » 1 (6 > a) it turns out. 


V 


A{b — a) 


52 


X 1^1 


b-2a-l 


(12a) 


(12b) 


and to guarantee ry <C 1 one needs a < (6 — l)/2. 

Therefore, each of the following inequalities guarantee 
adiabaticity in the very neighborhood of t = 0: 


where 


fk = ‘fk{0~,-oo ), 

‘fit = ¥5fe(+oo,0+). 

Now, consider what follows: 

(■*/'±(s)IV’±(s)) = 0, 

(■*/'±(s)|V’t(s)) = 


a> 2b+l, (13a) 

6 > 2a-hi. (13b) 

The last case to consider is given by the situation |a| ~ 
1, which deserves a specific analysis since it requires the 
analysis of further terms in the Taylor expansion of both 
A(t) and fl(t). 


and, more important, 

{^+(0+)|V'-h(0")) = COS0, 

(^_(0''‘)|'!/'-(0“)) = COS0, 

(V'-h(0+)|V'-(0")) = -sin0, 
(V^-(0 +)|'!/'-h( 0“)) = sin0. 


(19a) 

(19b) 

(19c) 


(20a) 

(20b) 


(21a) 

(21b) 

(21c) 

(21d) 
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with 

0 = 0 ( 0 +)- 0 ( 0 "). ( 22 ) 

From Eqs. (EH we can derive the probabilities of the 
nonadiabatic transitions. The failure of adiabatic pas¬ 
sage in the overall process can be quantified through the 
parameter 0: 

l(l/’±(+oo)| 17(-l-oo, —oo) \'tp±{—oo))\^ = cos^ 0 , (23a) 

|(^±(-|-oo)| 17(-|-oo, —oo) |'!/':f:(—oo))|^ = sin^ 0 . (23b) 

The first of such two equations, Eq. (I23al) . says that 
each of the two eigenstates of H, \tp+{t)) or fol¬ 

lows the adiabatic change of the Hamiltonian with prob¬ 
ability cos^ 0. On the other hand, Eq. (I23bl) says that 
each of such two states undergoes a transition to the other 
state with probability sin^ 0. 

The nonadiabatic transitions are determined by the 
abrupt change of the eigenstates of the Hamiltonian at 
t = 0. It is interesting to point out that, from the dynam¬ 
ical point of view, nothing really happens around t = 0, 
because the Hamiltonian is essentially zero in that time 
interval. Nevertheless, what happens from there on is 
mainly determined by the possible abrupt change of the 
Hamiltonian eigenstates occurring at t = 0. 

The case 0 = 0 can be obtained when a —>■ 0, which 
implies 0(0+) = 0(0“) = 7r/4. On the contrary, the case 
0 = tt/2 (—7r/2) can be obtained when a —>■ —oo (d-oo) 
on the left of t = 0 and a —-boo (—oo) on the right. 
Therefore, we can say that when A is infinitesimal of 
higher order than H there is no abrupt change of the 
eigenstates and then there is a total transfer of popula¬ 
tion. On the other hand, when H is infinitesimal of higher 
order than A we expect an abrupt change of the eigen¬ 
states and then a very low population transfer efficiency. 
When A and ft are infinitesimals of the same order, a 
specific analysis of the behavior of 0 is necessary. 

All these results have been derived under the assump¬ 
tion that the adiabatic approximation is valid everywhere 
far from t = 0, and that a possible critical behavior can 
occur only when the levels cross (t = 0). 

In the standard treatment of Landau-Zener- 
Stueckelberg-Majorana processes, fl is assumed to 
be constant and A is considered as a linear function of 
time. In our examples, both for the sake of simplicity 
and to make it easier to compare our situations with the 
standard one, we have always assumed A = nt, with n 
being a suitable parameter. 

In Eig. m we show the population transfer efficiency 
(i.e., the final population of the initially unpopulated 
state) when the H-pulse is essentially constant except 
for the neighborhood of the level crossing: H = Ho(I — 
exp(—|t/(T|)). The parameter a measures the size of the 
region where ft mainly varies. In Fig. [5] we consider a 
case where the shape of ft is far from being constant: 
ft = Aflo sin(7rt/T), where T represents the duration of 
the whole process, while A determines the maximum val¬ 
ues of the pulse. In both Figs. [T] and [21 a very good 



(a) 



(b) 


FIG. 1: (Color online) Frame (a) shows the asymptotic popu¬ 
lations of the bare state |t} (red solid line) and the theoretical 
prediction of the efficiency (blue dashed line), cos^ 0, vs the 
parameter a (in units of 2/flo). Here ft = f2o(l— exp(—|t/(T|)), 
A = Kt, with k/Q,Q = 0.025 and kT/Qq = 50, where 2T is the 
duration of the experiment. Frame (b) shows the coupling 
Q (red solid line) and the detuning A (blue dashed line) for 
a — 2/Ho. The amplitudes of the pulses are in units of Ho/2 
and time is in units of 2/Ho. 


agreement between the results from the numerical anal¬ 
ysis and the theoretical predictions (cos^ 0) is evident. 

Figure [3] compares numerical and analytical results for 
ft = Xf3tsva{'Kt/T). We observe nearly zero population 
transfer efficiency, which is due to the fact that a = A/ft 
diverges at t = 0, assuming negative values for t < 0 and 
positive values on the other side, which makes 0(O“) = 0 
and 0(0+) = 7r/2. The numerical results and theoretical 
predictions have a very good agreement also this time. 
Of course, in order to support the applicability of our 
theoretical analysis, we have always carefully checked 
(through numerical evaluation of ry) that the adiabatic 
approximation is valid essentially everywhere, with the 
only possible exception of t = 0 and its vicinity. 


B. Transitions around the level crossing 

It can happen that the validity of the adiabatic ap¬ 
proximation is broken before reaching the level crossing 
and for a while after it has occurred. In such a case, 
assuming that it is valid out of the (small) time interval 
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FIG. 2: (Color online) Frame (a) shows the asymptotic pop¬ 
ulations of the bare state |t) (red solid line) and the theo¬ 
retical prediction of the efficiency (bine dashed line), cos^ 0, 
vs the parameter A (dimensionless). Here Q = Aflo sin(7rt/r), 
where 2T is the duration of the experiment, and A = Kt, with 
k/Q,o = 0.025 and kT/Qq = 5. Frame (b) shows the conpling 
Q (red solid line) and the detuning A (blue dashed line) for 
A = 1. The amplitudes of the pulses are in units of flo/2 and 
time is in units of 2/flo. 


FIG. 3: (Color online) Frame (a) shows the asymptotic pop¬ 
ulations of the bare state |t) (red solid line) and the theo¬ 
retical prediction of the efficiency (blue dashed line), cos^ 0, 
vs the parameter A (dimensionless). The inset is a ‘zoom’ 
that shows the same quantities in a different y-range. Here 
fl = Xf5tsm{nt/T), where 2T is the duration of the experi¬ 
ment, and A = ret, with = 5 and /3T^ = 200. Frame 
(b) shows the coupling fl (red solid line) and the detuning A 
(blue dashed line) for A = 1. The amplitudes of the pulses 
are in units of \/p/2 and time is in units of \/2/p. 


[—T, r], one can write: 

[/(-boo, —oo) = [/(-boo, T)f/(r, —r)C/(— r, —oo) 

« C/A(+oo,T)17(r,-r)C/A(-T,-oo). 

(24) 


Since r is very small and assuming that H does not ex¬ 
hibit singularities at any time, we can make the following 
(rough) approximation: 


U{t,-t) 


I-iy H{s)ds, 
exp i J H{s)ds^ . 


(25) 


If A is an odd function, then one easily obtains: 


could be roughly given by: 

I (V’iMI C/(T,-r) \i^±i-T)) p 
~ cos^ ^ cos^ [ 61 (t) — d(—r)] 

-b sin^ ^ sin^ [6»(r) -b 6»(-r)]. (27) 

If the adiabaticity is broken far from the level cross¬ 
ing then all this treatment is not able to describe the 
behavior of the system. On the other hand, even under 
the proper hypotheses about adiabaticity, Eq. dUl) may 
give a wrong prediction because the first order approxi¬ 
mation to evaluate t/(—r, r) can be insufficient. In this 
case, further terms of the relevant Dyson expansion have 
to be used. 


U (r, —r) ~ cos p I— i sin ^ , 



D(s)ds. (26) 


Therefore, assuming adiabaticity everywhere out of the 
interval [—r, r], the efficiency of the population transfer 


V. DISCUSSION 

In this paper we have analyzed the dynamics of a two- 
state system whose energies undergo a level crossing. 
This is a kind of situation that cannot be treated by the 
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LZSM model, because the interaction between the two 
states approaches zero when the bare energies cross. 

The theoretical analysis makes it possible to predict 
the transition probability from one bare state to the other 
on the basis of the way the bare energies and the inten¬ 
sity of the interaction between them approach zero at 
the same time. Indeed, under the assumption that the 
adiabatic approximation is valid essentially everywhere 
except, possibly, around the crossing, the efficiency is 
mainly affected by the possible abrupt changes of the 
eigenstates of the Hamiltonian. In fact, if the bare ener¬ 
gies cross faster than the interaction there is no abrupt 
change of the eigenstates and the population transfer ef¬ 
ficiency tends to unity, while in the opposite case there is 
the most dramatic abrupt change of the eigenstates (each 
eigenstate immediately after the crossing is orthogonal to 
the corresponding eigenstate just before the crossing) and 
the efficiency tends to zero. If A and H approach zero as 
infinitesimal of the same order, then a specific analysis 
is required and the exact value of their ratio determines 
the efficiency of the population transfer. 

Our theoretical analysis is supported by numerical cal¬ 
culations. We have always found a good agreement be¬ 
tween the theoretically predicted efficiency and the corre¬ 


sponding numerical result, and in some cases this agree¬ 
ment is perfect. 

It is worth noting that in the very neighborhood of the 
level crossing nothing really happens, from the dynamical 
point of view, in the sense that the state of the quantum 
system essentially does not evolve (being H Ri 0), but 
what happens at the level crossing significantly affects the 
dynamics from there on. This circumstance is typical of 
any sudden change of the Hamiltonian, but in our case 
the peculiarity is that such a change is determined by 
infinitesimal quantities and the way they compare each 
other. Somehow, it resembles a symmetry breaking, since 
the perfectly symmetric situation at the crossing {H = 
0) is significantly altered, with important consequences 
on the system, by very small perturbations such as the 
pulses around t = 0. 
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